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ABSTRACT 

The present determination of the absolute magnitude My(RR) of RR Lyrae stars 
is twofold, relying upon Hipparcos proper motions and trigonometric parallaxes sepa- 
rately First, applying the statistical parallax method to the proper motions, we find 
< My(RR) >= 0.69 ± 0.10 for 99 halo RR Lyraes with <[Fe/H]> =-1.58. Second, 
applying the Lutz-Kelker correction to the RR Lyrae HIP95497 with the most accu- 
rately measured parallax, we obtain My(RR) = (0.58-0.68)lo;3i at [Fe/H]=-1.6. Fur- 
thermore, allowing full use of low accuracy and negative parallaxes as well for 125 RR 
Lyraes with - 2.49<[Fe/H]<0.07, the maximum likelihood estimation yields the rela- 
tion, M y (RR) = (0.59±0.37) + (0.20±0.63)([Fe/H]+1.60), which formally agrees with the 
recent preferred relation. The same estimation yields again < My(RR) > = 0.65 ±0.33 
for the 99 halo RR Lyraes. Although the formal errors in the latter three parallax es- 
timates are rather large, all of the four results suggest the fainter absolute magnitude, 
My(RR)^0.6-0.7 at [Fe/H]=-1.6. The present results still provide the lower limit on 
the age of the universe which is inconsistent with a flat, matter-dominated universe and 
current estimates of the Hubble constant. 

Subject headings: astrometry — stars: luminosities — stars: kinematics — stars: statis- 
tics 
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1. INTRODUCTION 

The absolute magnitude of RR Lyrae stars My (RR) 
is an indicator of distance and age in the Galaxy. In 
particular it is a key parameter for determining the 
age of Galactic globular clusters (GCs) which should 
place a lower limit on the age of the universe. How- 
ever, the values of My (RR) reported so far defies a 
consensus and splits into faint and bright values. 

The faint value of My (RR)=0.70-0.75 at the char- 
acteristic halo metallicity [Fe/H]=-1.6, which gives a 
shorter distance scale or larger age of clusters, is de- 
rived from halo RR Lyraes in the fields by the statis- 
tical parallax method (Hawley et al. 1986; Barnes & 
Hawley 1986; Strugnell, Reid, & Murray 1986; Lay- 
den et al. 1996) and by the Baade-Wesselink analysis 
(Liu & Janes 1990a; Jones et al. 1992). On the other 
hand, the bright value of M y (RR)=0.44 at [Fe/H]=- 
1.9, which gives a longer distance scale or smaller age, 
is derived from RR Lyraes in the GCs of the LMC, 
assuming a distance modulus to the LMC to be 18.5 
(Walker 1992). Sandage's pulsation theory (Sandage 
1993) and RR Lyraes in the metal-poor GCs of M5 
and M92 (Storm, Carney, & Latham 1994) also sup- 
port the bright value, but this is not the case for RR 
Lyares in M4 (Liu & Janes 1990b). 

Such a dichotomy in current estimates of My (RR) 
is the dominant source of uncertainty in the age de- 
termination of GCs. The faint My (RR) raises a well- 
known problem such that the age of GCs exceeds the 
age of the universe if the Hubble constant is adopted 
as H = 60 — 80km s _1 Mpc -1 , and this is the central 
argument in favor of nonzero cosmological constant 
(Bolte & Hogan 1995). However, a possibility of sway- 
ing back to the bright My (RR) has recently been dis- 
cussed since the preference of its bright value was in- 
ferred from the Hipparcos measurements of parallaxes 
of Cepheids (Feast & Catchpole 1997) and subdwarfs 
(Reid 1997). 

The astrometric observations of field RR Lyraes 
have been available to us as a part of the Hippar- 
cos programme assigned to the proposers in 1982. 
This paper is therefore the first report in the main 
journal on the direct determination of My (RR) from 
statistical treatments of Hipparcos proper motions 
and trigonometric parallaxes of RR Lyraes (cf., Tsuji- 
moto, Miyamoto, & Yoshii 1997; Fernley et al. 1997). 
The derived results from the proper motions in §2 
and from the parallaxes in §3 have confirmed the faint 
My (RR), thus suggesting that the so-called age dis- 



crepancy would not be resolved or alleviated yet. 

2. My (RR) DERIVED FROM HIPPARCOS 
PROPER MOTIONS 

The Hipparcos proper motion system is inertial, 
on a global sense, within the error ± 0.25 mas/yr and 
the mean error (cx^h of the Hipparcos proper motions 
[(/•Oh) (/As)h] for the RR Lyraes considered here is 
[±3 mas/yr, ±2 mas/yr]. Thus, the previous proper 
motions are remarkably improved by the Hipparcos 
mission. The most essential part of improvements of 
the previous proper motions [(/i^) p , {^s)p\ may be 
described by the rotation vector f2 of the previous 
proper motion system with respect to the Hipparcos 
system. The vector components (O x , £l y , fl z ) in the 
equatorial rectangular coordinates are derived from 
the equations of condition: 

[(/•Op _ (aOh] = -fix sin (5 cos a - fi,, sin 5 sin a 
[(Hs) P ~ (m<s)h] = +f2 x sina - n y cosa 



+n z cos s 




Applying eq.[l] to the proper motions of 99 halo 
RR Lyraes with <[Fe/H]>= -1.58, which are com- 
mon to Table 2 compiled by Layden et al. (1996) and 
Hipparcos Catalogue, we have 

n x = +1.61 ± 0.15mas/yr } 

fl v = +4.31 ± 0.13mas/yr I . (2) 

Q z = -1.27 ± 0.13mas/yr J 

Thus, the previous proper motions are systemat- 
ically improved by the amount of \ft\ ~ 5 mas/yr 
larger than the mean error (ct p )h in Hippracos proper 
motions for the RR Lyraes. Encouraged by the clear 
improvement indicated by eq. [2] , we apply the statis- 
tical parallax method to Hipparcos proper motions to 
improve My(RR) previously determined. 

The method is a direct application of well-understood 
kinematics of sample stars to the observations (proper 
motions, radial velocities, and apparent magnitudes), 
when a single luminosity class of sample stars can 
be specified within a small intrinsic magnitude dis- 
persion. The present maximum likelihood statistical 
analysis of RR Lyraes is similar to the one done on 
RR Lyraes by Hawley et al. (1986). Introducing a 
model that predicts the spatial motion of each star at 
each position in the Galaxy, we estimate the model 
parameters by the maximum likelihood method. The 
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Table 1: Results of statistical parallax analysis 



Altars <[Fc/H]> Uq Vq Wq a v a v a w My (RR) 



Hawley et al. (1986) 77 +21±19 

Layden et al. (1996) 162 -1.61 -9±14 

Fernley et al. (1997) 69 -1.66 

this work 99 -1.58 -12±17 



-184±17 -4±11 166±76 114±52 91±40 0.73±0.18 
-210±12 -12±8 168±13 102±8 97±7 0.71±0.12 

0.73±0.18 

-200±11 +2±3 161±13 105±9 88±7 0.69±0.10 



present model includes a reflex solar motion with re- 
spect to the galactic center and an ellipsoidal velocity 
distribution for RR Lyraes, but neglects the differ- 
ential galactic rotation. Assuming that the velocity 
residuals u, the observed velocity minus the velocity 
expected from the model, follow a Gaussian distribu- 
tion with zero mean, we have the logarithmic likeli- 
hood InL of the aggregate of all the residuals: 

lnL = -^T l \n{\M\ + u'M~ 1 v) + constant , (3) 

where M is the covariance tensor given by the expec- 
tation M= <uu'>. 

The model parameters describe the velocity com- 
ponents (Uq,Vq,Wq) of the reflex solar motion with 
U, V, and W axes pointing to the galactic cen- 
ter, the galactic rotation, and the north galactic 
pole, respectively, the velocity dispersion components 
{au,ov ,o~w), the covariances (ajjv ,o~uw -,o~vw), the 
absolute magnitude My (RR), and the intrinsic mag- 
nitude dispersion um in My (RR) . The observational 
errors uy r in the radial velocity and (a^ , cr^ s ) in the 
proper motions for each star are taken from Layden 
(1996) and Hipparcos Catalogue, respectively. The 
intrinsic dispersion om is parameterized by o^, the 
dispersion in the distance scale parameter k. Here we 
assume (7^=0.1 (Layden et al. 1996). The present 
model therefore incorporates 10 parameters, which 
are solved by the numerical minimization technique 
of simplex optimization. 

Combining Hipparcos proper motions with the 
data of [Fe/H], radial velocity, apparent magnitude, 
and interstellar extinction taken from Layden et al. (1996) 
and Layden (1994), we have carried out the present 
analysis for 99 halo RR Lyraes with < [Fe/H] > =-1.58, 
and found that all of the covariances {(Juv ,o~uw ,cvw) 
are nearly equal to zero. Therefore, we have re- 
peated the similar analysis to determine 7 parame- 



ters, putting all of the covariances equal to zero. The 
results are shown in Table 1. We find < My(RR) > 
= 0.69±0.10 together with the solar motion with re- 
spect to the galactic center (Uq,Vq,Wq) = (-12±17,- 
200±ll,+2±3) km s _1 and the velocity dispersions 
{au,a v ,a w ) = (161±13,105±9,88±7) km s" 1 . The 
RR Lyrae stars as a whole rotate with V rot ~30km 
s , assuming the local standard of rest Vlsr=220 
km s _1 and the solar motion of 16.5 km s _1 in the 
direction ^=53° and 6=25°. These results are in good 
agreement with those recently obtained by Hawley 
et al. (1986), Layden et al. (1996), and Fernley et 
al. (1997) (see Table 1) within the error estimates. 
Thus, the absolute magnitude of RR Lyraes derived 
from Hipparcos proper motions, < My >« 0.70 at 
<[Fe/H]> ?y -1.6, confirms the recent results from 
the ground-based observations. 

In general, a linear relation between My (RR) and 
[Fe/H] has been assumed. The slope of My (RR)- 
[Fe/H] relation has a strong influence on the inferred 
age differences among the globular clusters, which 
give a measure of the dynamical timescale for halo 
formation. The Baade-Wesselink results give the 
values of My (RR) over the range of [Fe/H]= -2.2 
to 0, which are fitted by the relation: My (RR) = 
(1.02±0.03) + (0.16±0.03)[Fe/H] (Jones et al. 1992) or 
My (RR) =(1.04±0.10) + (0.21±0.05)[Fe/H] (Skillen et 
al. 1992). However, the controversial situation holds 
also for determining the slope of My (RR)-[Fe/H] re- 
lation. Sandagc's pulsation theory gives much steeper 
slope, whereas statistical parallax analyses (Hawley 
et al. 1986; Strugncll et al. 1986; Layden et al. 1996) 
support very weak (or no) dependence of My(RR) 
on [Fe/H]. It is, of course, possible to incorporate the 
linear relation between My(RR) and [Fc/H] into the 
present analysis by dividing the available sample into 
bins of [Fe/H], but the available number (26) of Hip- 
parcos RR Lyraes with [Fe/H] > - 1.0 is too small 
to yield the reliable relation over the wide range of 
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[Fe/H] . 

3. My(RR) DERIVED FROM HIPPARCOS 
PARALLAXES 



Fig. 1. — The comparison between Hipparcos paral- 
laxes 7Th of 125 RR Lyraes and the photometric par- 
allaxes 7T p (see eq.[5]). The straight line indicates the 
relation 7TH=7r p . 

The relative error a n /irn for 173 Hipparcos RR 
Lyraes with the standard error a n of Hipparcos par- 
allaxes 7Th distributes such that each number of stars 
with ov/ttjj less than 1 and larger than 1 is 51 and 
65, respectively. The remaining 57 stars have even 
negative parallaxes. Thus, almost all of the par- 
allaxes are measured with very large error, except 
one case (HIP95497=RR Lyr) with ov/tt h = 0.135. 
In Figure 1 Hipparcos parallaxes 7Th are compared 
with the photometric ones 7r p defined by r p (kpc) = 
10 o.2(v-M v (RR)-io) ) where y ; s trie apparent magni- 
tude corrected for interstellar extinction and My(RR) 
the calibrated absolute magnitude given by the rela- 
tion 

My (RR) = a + 6([Fe/H] - c) (4) 

with a=0.91, 6=0.20, c=-1.60 (sec eq.[5]). The filled 
and open circles in the figure indicate the parallaxes 
with ov/7Th smaller than 1 and others, respectively. 
It is noticed that the individual deviation of Hippar- 
cos parallaxes from the photometric ones is extremely 



large, but Hipparcos parallaxes including low quality 
and negative ones as a whole agree with the photo- 
metric ones. 

The statistical properties of Hipparcos parallaxes 
of their errors have been investigated extensively by 
Arenou et al. (1995), who have shown that Hipparcos 
parallaxes are free from any global zero-point error up 
to ±0.1 mas and their errors are normally distributed 
with zero mean. This means that Hipparcos paral- 
laxes including even the negative ones for any ensem- 
ble of stars are statistically well-behaved. This is the 
case for the ensemble of RR Lyraes as well, as is shown 
in Figure 1. These statistical regularities of Hippar- 
cos parallaxes suggest a relevant technique for esti- 
mating statistically the absolute magnitude My(RR) 
of RR Lyraes to be applicable to the ensemble of RR 
Lyrae parallaxes, even though each individual paral- 
lax measured is largely erroneous, giving negative one 
as well. Here we attempt to derive My(RR) directly 
from Hipparcos parallaxes, relying upon two statisti- 
cal methods to find the unbiased estimate. 

a) Lutz-Kelkcr Correction 

Lutz & Kclkcr (1973) have shown that only when 
the relative error ov/V is smaller than 0.175, an un- 
biased estimate of the absolute magnitude can be re- 
trieved from the measured parallaxes 7r'. We have 
only one case of HIP95497 with ov/th = 0.135, to 
which Lutz-Kelkcr correction is applicable. The un- 
corrected absolute magnitude of this star is My (RR) 
= 0.82+ "_ 31 , which is derived formally from the mea- 
sured parallax. On the basis of Hanson's (1979) for- 
mulation to estimate the correction, we have the cal- 
ibrated absolute magnitude range My(RR) = (0.62- 
0.72)±°;^ for HIP95497 ([Fe/H] =-1.37), correspond- 
ing to the assumed spatial density of RR Lyraes 
n(r)=const. and n(r) oc r~ 2 as the lower and up- 
per bounds, respectively. If we assume 6=0.20 in 
eq.[4], the above range of My (RR) for HIP95497 cor- 
responds to My (RR) = (0.58-0.68)1°'^ at [Fe/H]=- 
1.6. 

b) Maximum Likelihood Calibration 

In order to retrieve statistically an unbiased es- 
timate of the absolute magnitude together with its 
intrinsic dispersion for a single luminosity class of 
magnitude limited stars from their measured paral- 
laxes with errors, Smith (1987,88) has proposed a 
maximum likelihood principle for correcting simulta- 
neously the Malmquist bias and the Lutz-Kelker bias, 
allowing full use of low accuracy and negative paral- 
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laxes. 

In line with the above principle we attempt to find 
the best esimate of the absolute magnitude My (RR) 
and of the intrinsic magnitude dispersion ctm for 
125 Hipparcos RR Lyraes in the metallicity range - 
2.49<[Fe/H]<0.07. We use here the numerical algo- 
rithm developed by Ratunatunga & Casertano (1991). 
Maximizing the likelihood function given by the prod- 
uct of the probabilities p(tth) of all the RR Lyraes 
considered, we obtain the best estimate of cm, a 
and b in eq.[4], where we adopt c = -1.60, taking 
into account the fact that the present distribution 
of RR Lyraes shows a peak at [Fe/H]«-1.6. Set- 
ting the allowable range of the measured parallaxes 
as [ 7r iowcr J 7r uppcr] = [^ 00 ,+ 00 ], an d assuming the num- 
ber density of RR Lyraes in space to be constant, we 
find the My (RR)- [Fe/H] relation: 

My (RR) = (0.59±0.37) + (0.20±0.63)([Fc/H] + 1.60) 

(5) 

with cr M = 2.6 x KT 4 ± 0.29. It is noticed that al- 
though the errors in the parameters estimated are 
large, the present My (RR)-[Fe/H] relation shows a 
formal agreement with the preferred relation My (RR) 
0.98 + 0.20[Fe/H] proposed by Chaboyer, Demarque, 
& Sarajedini (1996) without quoting error estimates. 
The insensitivity of the relation given by eq.[5] has 
been checked by putting um — 0.1-0.2 (Layden et 
al. 1996), [blower ,7r U pper]=[0,+oo], and n(r) oc r~ 2 
(Reid 1997), respectively. 

The large standard errors in eq. [5] are owing to the 
small numbers of degrees of freedom in the present 
statistics. To increase the freedom, putting 6=0 in 
eq.[4], we apply the same method to the 99 halo 
RR Lyraes with < [Fc/H] > =-1.58, whose absolute 
magnitude has been determined already on the ba- 
sis of the statistical parallax method. Then, we find 
< My (RR) >=0.65±0.33 and <7 M =4.1 x 10~ 6 ± 0.39 
at <[Fc/H]>= 1.58. This absolute magnitude can be 
compared with My (RR)=0.59±0.37 at [Fe/H]=-1.6 
given by eq.[5]. The validity of the present statistical 
analysis has been further checked by inspecting the 
distribution (7Th —TT p )/(J 7r for the 125 RR Lyrae stars. 
The median of the distribution is shifted only by 0.08 
from the Gaussian. 

4. CONCLUSION 

The present derivation of the absolute magnitude 
of RR Lyrae stars is twofold, relying upon Hipparcos 
proper motions and trigonometric parallaxes. First, 



applying the statistical parallax method to Hipparcos 
proper motions of 99 halo RR Lyraes, we have found 

< My (RR) >= 0.69±0.10 at < [Fc/H] > = -1.58 . 

(6) 

Second, applying the Lutz-Kelker correction to the 
most accurate parallax of HIP95497 (RR Lyr) with 
cr 7r /7rH=0.135 and [Fe/H]=-1.37, we have the magni- 
tude range M y (RR)=(0.62-0.72)^j^, which yields 

My(RR) = (0.58 - 0.68)^3? at [Fe/H] = -1.6 , 

(7) 

given the slope 0.20 of the My(RR)-[Fe/H] rela- 
tion. The lower and upper bounds in eq.[7] corre- 
spond to the assumed spatial density of RR Lyraes 
n(r)=const. and n(r) oc r~ 2 , respectively. Further- 
more, applyng the maximum likelihood algorithm 
to Hipparcos parallaxes including low accuracy and 
negative ones as well for 125 RR Lyraes with - 
2.49<[Fe/H]<0.07, we have found the relation (5) 
which gives 

My (RR) = 0.59 ±0.37 at [Fc/H] = -1.6 . (8) 

Applying the same method to the 99 halo RR Lyraes, 
we have found 



< My (RR) >= 0.65±0.33 at < [Fe/H] > = -1.58 . 

(9) 

Although the trigonometric parallax results have large 
errors, these results are consistent with the statisti- 
cal parallax result. All of the four results suggest 
the fainter absolute magnitude, My (RR)f«0.6-0.7 at 
[Fe/H]=-1.6. Note that three estimates from paral- 
laxes eqs.[7]-[9] are not completely independent, in 
regorous sense. 

Recently, Feast & Catchpole (1997) find My (RR) 
= 0.25 at [Fe/H] =-1.9, combining the LMC distance 
modulus of 18.7 determined by Hipparcos parallaxes 
of Cepheids with the data of LMC RR Lyraes by 
Walker (1992). On the other hand, Ried (1997) es- 
timates the distance scale to the globular clusters by 
subdwarf main-sequence fitting on the bais of Hip- 
parcos parallaxes, giving My(RR)=0.15 at [Fe/H]=- 
2.1. These bright values of My (RR) imply much 
younger age for the oldest globular clusters than the 
previous estimates (Bolte & Hogan 1995; Chaboyer et 
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al. 1996). However, the direct estimates of My (RR) 
in the present study support much fainter value than 
the indirect estimates obtained by them. The present 
results give My (RR) w 0.54 - 0.64 at [Fe/H]=- 
1.9, given the slope 0.20 of the My (RR)-[Fe/H] re- 
lation. Note that the Walker's My(RR) combined 
with the LMC distance modulus of 18.37 (Gould 
1995) derived from the SN 1987A ring distance yields 
My (RR)=0.57 at [Fe/H]=-1.9, which favors the present 
results. Our findings imply that the age of the old- 
est globular clusters in our Galaxy still conflicts with 
the standard cosmological model of a fiat, matter- 
dominated universe with the Hubble constant Hq = 
60 — 80kms _1 Mpc -1 estimated from almost all ob- 
servations (Bolt & Hogan 1995). It appears that the 
age problem remains still unresolved. 

This work has been supported in part by the grant- 
in- Aid for Scientific Research (08640336) and Center 
of Excellence (COE) research (07CE2002) of the Min- 
istry of Education, Science, and Culture in Japan. We 
would like to thank H. Saio for many fruitful discus- 
sions. 

REFERENCES 

Arenou, F., Lindcgrcn, L., Frocschle, M., Gomez, A. 
E., Turon, C, Perryman, M. A. C, & Wiclen, R. 
1995, A&A, 304, 52 

Barnes III, T. G., & Hawley, S. L. 1986, ApJ, 307, L9 

Bolte, M., & Hogan, C. J. 1995, Nature, 376, 399 

Chaboyer, B., Demarquc, K., & Sarajedini, A. 1996, 
ApJ, 459, 558 

Chaboyer, B., Demarque, K., Kcrnan, P. J., & 
Krauss, L. M. 1996, Science, 271, 957 

ESA, 1997, The Hipparcos and Tycho Catalogues, 
SP-1200, vol.1, chap.3 



Hawley, S. L., Jeffreys, W. H., & Barnes III, T. G. 
1986, ApJ, 302, 626 

Jones, R. V., Carney, B. W., Storm, J., & Latham, 
D. W. 1992, ApJ, 386, 646 

Layden, A. C. 1994, AJ, 108, 1016 

Layden, A. C, Hanson, R. B., Hawley, S. L., Klemola, 
A. R., Hanley, C. J. 1996, AJ, 112, 2110 

Liu, T., & Janes, K. A. 1990a, ApJ, 354, 273 

Liu, T., & Janes, K. A. 1990b, ApJ, 360, 561 

Lutz, T. E., & Kelkcr, D. H. 1973, PASP, 85, 573 

Ratnatunga, K. U., & Casertano, S. 1991, AJ, 101, 
1075 

Ried, I. N. 1997, AJ, 114, 161 
Sandage, A. 1993, AJ, 106, 703 

Skillcn, I., Fcrnley, J. A., Stobie, R. S., & Jameson, 
R. F. 1993, MNRAS, 265, 301 

Smith, H. 1988, A&A, 198, 365 

Smith, H. 1987, A&A, 171, 336 

Storm, J., Carney, B. W., & Latham, D. W. 1994, 
A&A, 290, 443 

Strugnell, P., Reid, N., & Murray, C. A. 1986, MN- 
RAS, 220, 413 

Tsujimoto, T., Miyamoto, M., Yoshii, Y. 1997, in 
HIPPARCOS Venice '97 Symposium, ed. M. A. C. 
Perryman (ESA) in press 

Walker, A. R. 1992, ApJ, 390, L81 



Feast, M. W., & Catchpole, R. W. 1997, MNRAS, 
286, LI 

Fernley, J., Barnes, T.G., Skillcn, I. et al. 1997, in 
HIPPARCOS Venice '97 Symposium, ed. M. A. C. 
Perryman (ESA) in press 

Gould, A. 1995, ApJ, 452, 189 



Hanson, R. B. 1979, MNRAS, 186, 875 This 2 _ column p re p rint was pre parcd with the AAS IATeX 

macros v4.0. 



6 



